MODELING  ANO  TRANSPORT  IN  SEMICONDUCTOR  KVICESCU) 
Arizona  state  univ  teak  d  k  ferrv  nov  *7 
ARO-21037.  9-EL  DAA029-04H(-R14« 


UNCLASSIFIED 


A 188  760 


1 


1  -ii  jh'.k' J*' A'A'cW* A*>ka >»‘ .fa' .W*A' .h1  Jh <  fa*  *«»  4«« JL' ju * ,fa*>b’ 


OTIC  FILE  COE* 


SECURITY  CLASSIFICATION  OF  This  PACE  (■»?,•«  D>1«  Entmrmd) 

F  REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 

ARO  21057.5-EL 


2.  GOVT  ACCESSION  NO. 


read  instructions 

_ BEFORE  COMPLETING  FORM 

~  RECIPIENT'S  CATALOG  NUMBER 


TITLE  (and  Subtltla ) 


Modeling  and  Transport  in  Semiconductor  Devices 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

1  Sept  1984  -  31  Aug  1987 
Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 


AU  THORf*J 


8.  CONTRACT  OR  GRANT  NUMBER^*; 


■■ 


D.K.  Ferry 


.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Arizona  State  University 
Tempe,  AZ  85287 


I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Army  Research  Office 

Research  Triangle  Park,  N.C.  27709 


DAAG29-84-K-0146 


10.  PROGRAM  element,  phoject,  task 
AREA  6  WORK  UNIT  NUMBERS 


13.  NUMBER  OF  PAGES 

81 


I*.  MONITORING  AGENCY  NAME  »  ADORESS<l/  dlttarant  trom  Controlling  Oil lea)  IS.  SECURITY  CLASS,  (ot  LI 1/a  import) 

Unclassified 

I  S>.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTIOH  STATEMENT  (ot  thla  Rm port) 


Approved  for  public  release; 
Distribution  unlimited. 


DTIC 

ELECTE 

DEC  1  5  1987 


'7.  DISTRIBUTION  STATEMENT  (ot  tha  abauact  antarad  in  Block  20,  It  dlttarant  from  Report) 


18.  supplementary  notes 


Ihe  view,  opinions  and/or  findings  contained  in  the  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 

'9-  KEY  WORDS  (Continue  on  rarataa  aidm  It  nacaaaary  and  Identity  by  block  numbar) 

Laser  excitation,  semiconductors,  hot  electrons. 


ABSTRACT  ("Cocrttmj*  ea  nr*r«*  aide  ff  rrmcaeearj  mod  idarutty  by  block  n  am  bar) 

This  report  covers  studies  of  picosecond  laser  excitation  of  dense  electron-holt 
planners  in  semiconductors.  Non-equilibrium  phonons  and  the  electron  hole 
interaction  are  found  to  be  import  processes  in  the  equilibration  of  the  plasma. 


F  ORM  > 

*  JAM  73  *4/3  EDITION  OF  I  MOV  65  IS  ULiOLE  l  C 


SECURITY  CLA5SIFIC  ATION  OF  THIS  PAGE  (Wh~m  Dmtm  tnlarad) 


t  4^4 


ABSTRACT 
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-V- 

The  advent  of  time-resolved,  subpicosecond  laser  pulse  techniques 

•  ;  <  •  i  ■ 

has  allowed  us  to  experimentally  study  the  ultrafast  relaxation 

v 

process  in  the  far-from  equilibrium  electron-hole  plasma  in  GaAs . 
Several  theoretical  studies  have  been  carried  out  to  explain  these 
experimental  measurements.  The  novel  ensemble  Monte  Carlo  simulation 
technique  was  used  to  investigate  the  cooling  process  of  the  hot 
laser-excited  carriers,  considering  the  roles  of  the  carrier-phonon 
interaction,  electron-hole  interaction,  and  the  hot  phonon  effects. 
The  time-evolving  distribution  function,  which  is  built  self- 
consistently  in  the  ensemble  Monte  Carlo  method,  has  rendered  it 
unnecessary  to  postulate  the  shape  of  the  distribution  function  or  to 
assume  equal  temperatures  for  the  electrons  and  holes,  ^the  effects  of 
different  excitation  energies  and  excitation  levels  were  examined 
without  applying  external  electric  fields.  Screening  effects  were 
investigated  separately  and  electron-electron  and  hole-hole 
interactions  were  ignored  in  the  calculations  to  emphasize  the 
electron-hole  interaction  and  hot  phonon  effects.  Degeneracy  was 
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taken  into  account  while  the  excitation  level  was  greater  thaip-lxtu 

_3 

cm  .  The  study  showed  that  the  hot  phonon  effects  do  slow  down  the 
cooling  rates  as  expected  and  the  electron-hole  interaction  slows  the 
cooling  rates  at  low  excitation  levels  while  enhancing  it  at  high 


excitation  levels. 
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CHAPTER  1 


INTRODUCTION 


1.1  Objective 


The  investigation  of  the  ultrafast  dynamics  of  hot  photoexcited 
carriers  in  semiconductors  is  crucial  for  understanding  the  basic 
physics  of  scattering  processes  in  semiconductors  and  has  motivated 
many  recent  experimental  optical  studies.  The  processes  that  govern 
the  carrier  dynamics  are  best  investigated  by  first  generating  these 
particles  in  a  nonequilibrium  state  with  a  short  laser  pulse,  and  then 
observing  their  subsequent  relaxation  phenomenon.  The  study  of  the 
transient  dynamics  of  carriers  excited  by  subpicosecond  laser  pulses 
in  GaAs  has  been  investigated  extensively  in  the  last  fev  years. 
Improvement  in  time  resolution  and  sensitivity  of  the  subpi cosecond 
measurement  system  have  allowed  the  investigation  of  processes  such  as 
extraction  of  relaxation  times  corresponding  to  different  scattering 
processes  [1],  velocity  overshoot  of  semiconductors  [2,3],  cooling 
rates  of  hot  photoexcited  carriers  [A, 5],  lifetime  of  LO  phonons  [6- 
8],  screening  of  electron-phonon  interaction  [9]. 

Several  authors  have  carried  out  theoretical  studies  attempting 
to  understand  these  measurements  [11,15].  Vithout  modern  computers 
one  had  to  assume  a  rather  simple  distribution  function  (e.g.  drifted 
Maxwellian)  and  a  simplified  band  structure  in  order  to  make  the 
transport  calculations  tractable.  More  recently,  fast  computers  and 
the  development  of  the  efficient  ensemble  Monte  Carlo  techniques  have 


allowed  one  to  avoid  these  assumptions  and  use  more  realistic  band 
structures  without  constraints  on  the  distribution  function.  The 


objective  of  this  thesis  was  to  investigate  an  ensemble  Monte  Carlo 
technique  which  simulates  the  ultrafast  relaxation  processes  that 
occur  on  the  very  short  subpicosecond  time  regime  and  to  evaluate  the 
roles  of  electron-hole  interaction  and  hot  phonon  effects  on  the 
cooling  process  of  the  hot  photoexcited  carriers  in  semiconductors. 

1.2  The  Approach 

The  important  scattering  mechanisms  in  semiconductors  have  been 
identified  and  the  theory  of  these  processes  have  been  well  developed 
in  the  last  few  years.  The  ensemble  Monte  Carlo  approach  developed 
here  has  taken  account  of  the  electron-phonon,  hole-phonon,  and 
electron-hole  scatterings,  and  included  the  hot  phonon  effects  in  the 
simulation,  which  means  the  phonon  system  now  is  considered  to  be  in  a 
nonequilibrium  state.  The  band  structure  of  GaAs  considered  here 
consists  of  three  nonparabolic  valleys  (T.L.X)  in  the  conduction  band 
and  two  parobolic  heavy-hole  and  light-hole  bands  in  the  valence  band. 
Since  carrier-carrier  scattering  relies  strongly  on  the  carrier  den¬ 
sities,  the  scattering  mechnisms  involving  the  weakly-populated  light 
holes  were  ignored  in  this  calculation.  However,  this  does  not  mean 
that  the  light  holes  do  not  play  any  role  in  the  whole  cooling 
process.  The  investigation  of  the  impact  of  the  light  holes  needs 
further  consideration.  Since  ve  assume  the  carriers  being  simulated 


are  all  generated  by  the  laser  pulse  rather  than  from  the  external 


3 


doping,  the  elastic  impurity  scattering  was  ignored  in  this  calcula¬ 
tion,  while  the  electron-electron  and  hole-hole  scatterings  also  were 
ignored  up  to  this  stage,  because  of  the  necessity  of  modifying  the 
program  due  to  the  inclusion  of  the  hot  phonon  effect.  This  simula¬ 
tion  has  been  carried  out  for  several  different  excitation  levels  and 
excitation  energies. 


1.3  Literature  Review 


One  of  the  most  interesting  phenomena  of  the  hot  photoexcited 
relaxation  experiments  has  been  the  unexpectedly  slow  cooling  of  the 
hot  distribution.  It  has  been  observed  that  the  cooling  slows  as  the 
carrier  density  increases  [13).  The  first  experiment  investigating 
the  energy  transfer  in  an  optically  excited  electron-hole  plasma  was 
done  by  HOpfel,  Shah,  and  Gossard  (1986)  [14].  The  total  energy-loss 
rate  for  electron  is  found  to  be  much  larger  in  the  presence  of  a  hole 
plasma  due  to  the  energy  flow  channelled  by  the  electron-hole  interac¬ 
tion  which  transfers  energy  from  the  hot  electrons  to  the  cold  holes. 
The  holes,  after  optical  excitation,  never  obtain  a  temperature  as 
high  as  the  electrons  because  of  the  much  heavier  effective  mass  of 
hole  and  the  realistic  band  structure  of  the  material.  Asche  and 
Sarbei  [15]  estimated  the  electron-hole  interaction  and  showed  that 
the  more  convincing  assumption  of  different  carrier  temperatures  leads 
to  a  much  slower  cooling  rate  than  obtained  for  equal  carrier  tempera¬ 
tures  (16J.  Although  the  assumption  of  equal  carrier  temperatures  in 
Ptttz's  approach  is  not  very  reasonable  from  a  physical  point  of  view, 
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CHAPTER  2 


THEORETICAL  BACKGROUND 

2.1  Band  Structure  of  GaAs 

The  band  structure  of  a  crystalline  solid,  i.e.  the  energy- 
momentum  (E-k)  relationship,  is  generally  obtained  by  solving  the 
Schrodinger  equation  of  an  approximate  one-electron  problem.  Quite 
often,  the  band  structure  is  very  complicated  and  can  only  be  calcu¬ 
lated  semi-empirically  utilizing  a  computer.  A  variety  of  numerical 
methods  have  been  used  to  theoretically  calculated  the  energy  bands  of 
solids.  The  methods  used  most  often  for  semiconductors  are  the  or- 
thogonalized  plane-vave  method  [10],  the  pseudopotential  method  [26], 
and  ic  •  p  method  [27].  Figure  2.1  shows  the  calculated  band  structure 
of  GaAs.  including  the  spin-orbit  interaction  [28].  Fig.  2.2  il¬ 
lustrates  the  important  region  of  the  GaAs  band  structure  in 
transport.  The  top  of  the  valence  band  at  k=0,  the  T  point,  is 
degenerate  with  the  upper  valence  band,  termed  the  heavy-hole  banc, 
the  lover  valence  band,  termed  the  light-hole  band.  About  0.34  ev 
[29]  below  the  T  point  lies  the  spin-orbit  split-off  banc.  The  mini¬ 
mum  of  the  conduction  band  also  lies  at  k=0,  which  makes  GaAs  a 
direct-gap  semiconductor  material,  in  contrast  to  Si.  The  energy  gap 

E  of  GaAs  is  temperature-dependent  and  considered  to  be  1.51  ev  for 
£ 

T  <7 7 0 K ,  1.439  ev  for  T>77°K  here. 

As  shown  in  Fig.  2.2.  the  valence  band  can  be  approximately 
fitted  by  two  parabolic  bands  with  different  curvatures:  a  wider  band 
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with  smaller  —  (the  heavy-hole  band)  and  the  narrower  band  with 
3  k 


32E 


larger  —r—  (the  light-hole  band).  The  effective  mass  of  the  carriers 
3  k 


•w  n 

is  generally  a  tensor  quantity  with  components  m.^  defined  as  m.^ 


( wr1- 
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reality,  these  two  bands  are  warped  and  nonparabolic 


as  a  result  of  their  degeneracy  and  the  interaction  among  the  bands. 
The  reason  we  treat  the  valence  band  as  parabolic  is  that  at  the 
present  time  one  can  get  analytical  expressions  for  the  carrier- 
carrier  scatterings  only  for  parabolic  bands.  There  are  two  sets  of 
satellite  conduction-band  valleys  in  GaAs.  The  L-valleys  have  minima 
at  the  boundary  of  the  Briilouin  zone  in  the  (111)  directions  at  an 
energy  about  0.30  ev  [30,31]  above  the  central  valley.  The  X-valleys 
have  minima  at  zone  edges  in  the  (100)  directions  at  an  energy  about 
0.50  ev  [32]  above  the  central  valley.  Due  to  the  tunability  of 
lasers,  we  can  control  the  excitation  energy  to  determine  whether  or 
not  excitation  of  the  electrons  up  to  the  satellite  valleys  is 
occurring.  These  three  valleys  are  nonparabolic  and  anisotropic,  each 
requiring  a  longitudinal  and  transverse  mass  component  to  describe  it. 

The  E-k  relationship  for  a  nonparaboli^  band  takes  the  form  [12] 


h"k" 


E( 1  *  c£ ) , 


and 
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where  a  is  the  nonparabolici ty  factor.  For  a  parabolic  warped  band, 
the  E-k  relationship  becomes  [37] 


E(k) 


h2k2 

2m 


s(  40 , 


where 


g(e,  ♦)  =  A  ±  [  B2  +  ^-fk2k2+k2k2+k2k2l  1 
L  y  x  z  y  z)  J 


1/2 


The  average  effective  mass  for  such  a  band  is  given  by 


where  the  positive  and  negative  signs  are  for  the  light-hole  and 
heavy-hole  respectively.  For  GaAs,  A=-6.98,  B=-4.5,  and  C=6.2  [33]. 


2.2  Scattering  Mechanisms 


The  scattering  processes  that  are  included  in  calculating  the 
time-evolved  carrier  distributions  on  a  subpicosecond  time  scale  can 
be  classified  into  two  groups  [34]:  processes  in  which  the  scattering 
rates  are  independent  of  carrier  concentrations,  such  as  unscreened 
carrier-phonon  scattering,  and  processes  in  which  the  scattering  rates 
depend  strongly  on  carrier  densities,  such  as  carrier-carrier 
scattering.  In  general,  the  carrier-density-dependent  processes 


10 


consist  of  the  screening  effect  of  the  carrier-phonon  scattering  and 
the  carrier-carrier  scattering,  while  the  carrier-density-independent 
processes  including  the  acoustic  phonon  scattering,  polar  optical 
phonon  scattering,  non-equivalent  intervalley  phonon  scattering,  and 
equivalent  intervalley  phonon  scattering.  The  carrier-carrier  scat¬ 
tering  can  be  classified  as  inter-  and  intra-valley  electron-electron 
scatterings,  inter-  and  intra-band  hole-hole  scatterings,  and 
electron-hole  scattering.  Since  the  Coulomb  cross  section  decreases 
rapidly  as  momentum  transfer  increases,  the  scatterings  caused  by 
ioni2ed  impurities  or  other  carriers  are  generally  considered  only  for 
intravalley  or  intraband  transitions  [35]. 

According  to  time-dependent  first-order  perturbation  theory,  the 
transition  probability  per  unit  time  for  scattering  from  a  state  £  to 
a  state  due  to  a  perturbing  hamiltonian  H'  is  given  by  the  Fermi- 
Golden  rule 


)  =  jpi<itiH'  |itf>i  6(e-e'  >, 

where  E  and  E'  are  the  eigenenergies  of  the  initial  and  final  states 
respectively.  It  is  convenient  to  use  Fourier  transformation  to 
interpret  H',  which  gives  the  squared  matrix  element  the  following 
form  [36,37]: 

|  <K  |  H '  [  k '  >  |  2  =  V(q)G(f*.k'). 
thus  the  transition  probability  is 


where  V(q)  contains  the  dependence  on  Ak  of  the  squared  Fourier  trans¬ 
formation  of  the  interaction  potential  and  G(k,k')  represents  the 
overlap  integral 

G(it,ic' )  =  I  J  nu*  (r)u„(r)d?  |2. 
ce11  k'  k 

The  overlap  factor  equals  unity  for  parabolic  bands ,  while  it  is  less 
than  unity  if  nonparabolici ty  is  significant.  The  overlap  effect  was 
ignored  for  electrons  because  G(k,k')  in  this  case  is  very  close  to 
unity.  The  overlap  factor  was  taken  to  be  0.5  for  holes  [52]  in  the 
calculations.  Actually)  for  the  intraband  transition  of  holes,  the 
overlap  factor  is  taken  as  [48] 

G(it,it')  =  i(l*3cos2e), 

and  for  interband  transitions, 

G(it,it')  =  jsin2©, 

where  ©  is  the  angle  between  k  and  k' .  The  various  scattering 
processes  considered  in  this  calculation  will  be  discussed  in¬ 
dividually  in  the  following  two  sections.  The  scattering 
probabilities  per  unit  time  with  respect  to  the  carrier's  energy  are 
well  formulated  for  each  scattering  mechanism. 


2 


2.2.1  Formulations  of  Electron  Scattering  Rates 

Acoustic  Phonon  Scattering 

This  scattering  process  involves  the  electrons  being  scattered  by 
intravalley  acoustic  phonons  through  deformation  potential  interaction 
which  is  due  directly  to  the  strain  caused  by  the  acoustic  waves  in 
the  crystals.  The  energy  involved  is  small,  thus  we  can  treat  it  as 
being  approximately  elastic  in  nature.  The  squared  matrix  element 
[37]  is  given  as 
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where  2^  is  the  acoustic  deformation  potential,  p  the  crystal  density, 
q  the  phonon  wave-vector,  u  the  speed  of  sound,  N  the  Bose-Eins tein 
distribution  function,  and  the  ±  signs  represent  the  absorption  and 
emission  respectively.  Since  we  are  not  interested  in  a  lattice 

C 

temperature  above  300  K,  scattering  processes  involving  more  than  one 
phonon  are  neglected.  According  to  the  elastic,  equiparti tion-energy 
approximation,  the  scattering  rate  for  nonparabolic  bands  is  given  as 
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The  angles  of  the  final  wave-vector  can  be  determined  by  random  num¬ 
bers  r^  and  r2  because  this  process  completely  randomizes  the 
momentum,  and 


cos©  =  l-2r^, 
4>  =  2rtr2> 


Polar  Optical  Phonon  Scattering 


(2.1) 

(2.2) 


In  polar  semiconductors  such  as  GaAs,  the  polarization  of  the 
crystal,  arising  from  the  vibrations  of  the  longitudinal  optical 
branch,  results  in  an  additional  interaction  with  the  carriers  termed 
the  polar  optical  interaction.  The  energy  of  polar  optical  phonons  in 


intravalley  transitions  can  be  assumed  constant,  given  by  hw0  =  k.gT  , 
where  T  is  the  equivalent  temperature  of  the  phonons.  The  polar 
optical  scattering  due  to  TO  phonons  is  negligible  in  contrast  to  LO 


phonons.  The  squared  matrix  element  is  [37] 
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where  w0  is  the  longitudinal  optical  mode  frequency,  and  e  ,  ec  the 
dielectric  constants  for  infinite  and  zero  frequencies  respectively. 
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N0  becomes  q-independent  and  has  the  form 


where  T  is  the  lattice  temperature.  The  scattering  probability  per 
unit  time  is  given  as  [12] 
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where  E'  =  E  *  hw0,  y(E)  =  E(l+o£),  and 


P0(E,E') 
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The  total  scattering  rate  for  polar  optical  scattering  decreases  as 
energy  increases  as  shown  in  Fig.  2.3.  This  scattering  is  strongly 
anisotropic  and  the  angular  dependence  of  the  scattering  is  given  by 
[36] 


P(9)dG  = 


_ sin  6d  6 _ 

r(E)  +  y(E'  )-2/y(E)y(E'  ) cos 6 


Using  random  numbers  r  and  r.  uniformly  distributed  between  0  and  1, 
we  get 


■ 

cose  =  (Uf)-(U2f)r  /  f , 


E '  =  E  +  nw.  , 
1 


where  Z  is  the  number  of  equivalent  final  valleys,  and  N\  is  the 
population  of  the  phonons  involved  in  the  transitions.  All  states  in 
the  energy-conserving  sphere  are  equally  probable  to  be  the  final 
state  after  scattering. 


Non-eouivalent  Intervalley  Phonon  Scattering 


This  scattering  mechanism  is  responsible  for  negative  differen¬ 
tial  mobility  as  the  origin  of  the  velocity  overshoot  phenomenon 
observed  in  GaAs.  Electrons  transfer  between  the  central  valley  and 
one  of  the  outer  satellite  valleys  either  by  absorption  or  emission  of 
a  phonon.  The  squared  matrix  element  is  given  as  [37] 


D:^n  f.  1  n 
°  ^  2  choc.  .  v  q  2  “  2  J  ’ 


where  D.^  is  the  deformation  potential  for  scattering  from  ith  valley 
to  jth  valley  induced  by  a  phonon  of  energy  hoc.  ^  .  The  scattering  rate 
is  given  as  [35] 


P  .(E)  =  (2-1 : 
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m*3/:D2. 
1 1 
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where  a.  is  the  nonparabolici ty  iactor  of  the  jth  valley,  and  Af... .  is 


the  energy  difference  between  the  bottoms  of  the  initial  and  final 
valleys.  Since  the  intervalley  scattering  randomizes  the  momentum, 
there  is  no  preferred  direction  for  the  final  state.  These  rates  are 
shown  in  Figs.  2.4  and  2.5. 


Electron-Hole  Scattering 


At  high  carrier  densities,  the  collisions  of  the  carriers  with 
each  other  need  to  be  considered.  The  electrons  and  holes  are  assumed 
to  interact  among  themselves  via  a  screened  Coulomb  potential 
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V(r)  =  4iiTXP(-6r)’ 


where  e  is  the  dielectric  constant  of  the  material  and  3  is  the  in¬ 
verse  screening  length  which  can  be  calculated  self-consistently  in 
the  ensemble  Monte  Carlo  method.  According  to  the  Fermi-Golden  rule, 
the  transition  probability  for  two  carriers  to  transit  from  state 

*4  “4  -4  < 

(kpk^)  to  state  ( k  £ .  k  A )  due  to  the  coulomb  interaction  is 


where  f  and  E  are  the  occupation  probabilitv  anc  the  energv  at  state 
k  k 

The  matrix  element  is  giver. 
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vhere  the  initial  and  final  states  are  considered  to  be  plane  waves. 


The  scattering  rate  of  the  carrier-carrier  scattering  is  obtained  by 


*4  *4  ^4 

summing  over  k^ ,  k£  and  k^ .  Since  the  final  state 


degeneracy  is  considered  separately  in  the  EMC  method  151],  all  the 


final  states  are  assumed  to  be  unoccupied,  so  that,  we  can  set  f 


f^  =  0,  and  f  =1.  Under  these  conditions,  the  scattering  rate  for 


an  electron  in  state  |k^>  scattered  by  a  hole  in  state  |k^>  is  given 


as  [36] 
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After  special  manipulations  (for  the  details,  see  Ref. [36]),  Pg^(k^) 


can  be  expressed  as 
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(2.3) 


vhere  the  relative  wave  vector  g  and  reduced  mass  u  are  defined  as 
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Rewriting  (2.3)  in  integral  form,  we  obtain 
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where  Q  ,  =  2u I  —  -  —  I,  and  p  is  the  concentration 
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Equation  (2.4)  is  now  used  by  replacing  the  distribution 
integral  by  a  particle  sum  in  the  Monte  Carlo  procedure. 


2.2.2  Formulations  of  Hole  Scattering  Rates 

A  parabolic  and  spherically  symmetric  two-banc  mode 
sidered  when  calculating  the  scattering  rates  for  holes, 
split-off  band  is  omitted.  In  the  valence  banc,  the  polar  op 
deformation  potential  scatterings  dominate  while  the  impurity 
ing  is  neglected  in  this  calculation. 


Acoustic  Phonon  Scattering 


The  acous  t  i  c -phonon  coupling  constant  ~  :r  giver,  as  [fl  ] 
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where  s^,  s  are  the  longitudinal  and  transverse  sound  velocities,  a, 
b  and  d  are  the  deformation  potential  constants  [52].  The  average 
longitudinal  and  transverse  elastic  constants  c^  and  c£  are  given  as 

C1  =  5  (  3c11  +  2c12+4c44  )’ 

Ct  =  5  (  Cll_C12+3c44  )’ 

where  c^,  c^  and  c^  are  the  crystal  elastic  constants.  The  values 
of  all  the  parameters  used  in  the  calculations  are  shown  in  Table  2.1 
[53].  The  resulting  scattering  rate  for  intraband  transitions  is 
given  as  [52] 
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For  interband  transitions,  ve  have 


Bulk  material  parameters 

2 

Crystal  elastic  constants  (dyn/cm  ): 

11.88x10 

C12 

5.38x10 

C44 

5.49x10 

[100]  longitudinal  sound  velocity  (cm/sec) 

4.73x10 

[100]  transverse  sound  velocity  S  (cm/sec) 

3.34x10 

Scattering-rate  parameters 

Effective  masses:  heavy-hole  band  m^ 

0.45m 

o 

light-hole  band 

0.082m 

0 

Optical-phonon  energy  (eV) 

0.036 

Deformation-potential  constants  (eV):  a 

6.7 

b  1.7 

d  4.4 

Table  2.1  Parameters  for  GaAs  hole- transport  program 
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where  result  from  the  integration  of  G(0)  over  all  scattering 

angles  6,  and  it  is  found  G^  =  G^2  =  G^  =  G^^  =  0.5. 

Non-Polar  Optical  Phonon  Scattering 


The  scattering  rate  for  intraband  processes  is  given  by  [ 36 ] 
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while  for  interband  transitions,  it  becomes 
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The  angles  of  the  final  wave-vector  after  scattering  can  be  chosen  by 
(2.1)  and  (2.2)  because  this  scattering  process  completely  randomizes 
the  momentum.  These  rates  are  shown  in  Figs.  2.6  and  2.7. 

Polar  Optical  Phonon  Scattering 

The  scattering  rate  for  intraband  transitions  is  [36] 
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The  angular  dependence  for  intraband  transitions  is 


Pii(e)d6 


( l*3cos^  6)si  n  6d  9 
V+k'*  -2kk'cos9  ’ 


while  for  interband  transitions,  it  becomes 
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These  scattering  rates  are  plotted  in  Figs.  2.8  and  2 


Hole-Electron  Scattering 


The  scattering  rate  for 


an  hole  in  state  ! k. 


electron  in  state 


I  k„> 


is  given 


as  [36] 
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and  n  is  the  concentration  of  electrons. 

2.3  Ultrafast  Relaxation  Mechanisms 

Progress  in  subpicosecond  laser  pulse  techniques  has  allowed 
direct  probing  of  the  relaxation  processes  of  the  hot  photoexcited 
carriers  in  semiconductors.  These  Vorks  have  shown  that,  on  the 
picosecond  time  scale,  a  quasi -equilibrium  electron-hole  plasma  is 
established  through  the  combined  action  of  carrier-phonon  and  carrier- 
carrier  interactions.  Typically,  the  relaxation  of  hot  laser-excited 
carriers  can  be  divided  into  three  stages.  Primarily  mediated  by 
Frohlich  interactions,  there  is  a  very  fast  transfer  of  excess  energy 
from  the  carriers  to  the  LO-phonons  until  their  mutual  thermalization 
is  reached.  During  this  initial  stage,  the  amount  of  energy  received 
by  the  acoustic  phonon  system  is  negligible.  In  the  second  stage,  the 
hot  distribution  cools  towards  the  lattice  temperature  until  optical 
phonons  are  no  longer  active  in  removing  energy.  The  hybrid  system  of 
carriers  and  LO-phonor.s  releases  their  excess  energy  to  the  lattice  at 
a  rate  determined  by  the  process  of  LO-phonons  relaxing  their  own 
energy  to  the  lattice  through  the  anharmor.ic  effects,  while  the  car¬ 
riers  transfer  their  excess  energy  to  LO-phonons  to  maintain  a  mutual 
quasi-equilibrium.  The  characteristic  cooling  times  are  much  longer 
than  those  in  the  initial  stage.  The  last  stage  involves  the  final 
cooling  of  the  distribution  to  the  lattice  temperature  through  acous¬ 
tic  phonon  emission.  Sampaio  and  Luzzi  [38]  established  that  heating 


of  LO-phonons,  due  to  scattering  with  carriers,  begins  within  1  psec. 


the  energy  transfer  from  carriers  to  LO-phonons  occurs  over  some  5 
psec,  and  the  energy  transfer  from  carriers  to  acoustic  phonons  occurs 


on  a  nanosecond  time  scale.  Typically,  the  first  two  stages  in  the 
relaxation  last  more  than  100  psec.  While  the  subsequent  cooling  of 
this  electron-hole  plasma  towards  the  lattice  temperature  is  well 
studied  [39],  little  is  known  about  the  primary  stage  of  the  carrier 
relaxation,  and  the  role  of  electron-hole  scattering.  Recently,  there 
have  been  several  experimental  investigations  of  the  initial  stage  of 
the  ultrafast  relaxation  done  by  some  researchers,  such  as  Tang  and 
Erskine  et_  al_.  [40-42,34],  Oudar  and  covorkers  [43,44]  etc. 

In  the  initial  stage  of  the  excitation  process,  most  of  the 
kinetic  energy  is  given  to  the  electron  because  it  is  much  lighter 
than  the  holes.  Most  transitions  will  be  from  the  heavy-hole  band 
because  the  light  holes  have  a  much  lower  population.  Thus,  we  con¬ 
centrate  on  understanding  the  relaxation  of  the  electrons  in  the 
initial  stage.  The  electrons  cascade  down  through  the  conduction  band 
as  fast  as  they  can  until  their  mutual  thermalization.  Several  scat¬ 
tering  processes  will  play  important  roles  depending  or,  the  electron 
density.  At  very  low  densities  the  electrons  can  emit  LG  phonons 
successively  and  hop  down  the  bands  in  units  of  L0  phonon  energy.  At 
higher  electron  densities,  the  L0  phonons  and  plasmons  will  couple 
together,  the  hopping  now  is  in  units  of  the  coupled  mode  energy.  At 
very  high  electron  densities,  the  L0  phonon  interaction  will  be  fully 
screened  and  the  electrons  will  hop  down  the  conduction  band  in  units 
of  the  plasm, on  energy  [-1].  «e  should  note  that  scattering  of  the 


electrons  by  the  free  holes  is  less  efficient  than  the  electron- 


electron  scattering  in  relaxing  the  excess  energy  in  this  stage, 
because  most  holes  are  much  heavier  than  the  electrons.  So,  in  the 
initial  stage,  electron-electron  scattering  and  polar  optical  phonon 
emission  are  dominant. 

In  the  second  stage  of  the  relaxation,  most  carriers  have  inter¬ 
changed  most  of  their  initial  excess  energy  and  formed  a  quasi¬ 
equilibrium  distribution  around  the  band  edges.  As  discussed  above, 
hot  carriers  relax  basically  through  carrier-phonon  scattering  and 
carrier-carrier  interaction.  It  is  hard  to  elucidate  which  scattering 
occurs  more  frequently.  The  energy-loss-rate  is  a  sum  of  three 
processes:  LO  phonon  emission  by  both  electrons  and  holes,  and  TO 
phonon  emission  by  holes,  because  both  electrons  and  holes  can  con¬ 
tribute  to  the  phonon  emission  in  the  second  stage  [45].  The  carrier- 
acoustic  phonon  interaction  gives  a  minor  contribution  to  the 
relaxation  for  most  of  the  processes. 

2.2.1  The  Roles  of  Hot  Phonon  Effects  and  Screened  carrier-phonon 
interactions  in  relaxation 

Due  solely  to  these  carrier-phonon  scattering  mechanisms,  the  so- 
callec  "hot  phonon  effects"  may  arise,  as  observed  experimentally 
[46.13].  Vhen  excess  nonequilibrium.  LO  phonons  are  rapidly  generated 
in  the  initial  stage,  there  results  a  slow  carrier  relaxation  towards 
the  final  equilibrium,  with  the  lattice.  since  the  nonequi 1 i b r i urn 
phonons  are  re-absorbec.  As  the  generated  LO  phonons  are  larger  than 


the  equilibrium  phonon  population,  they  follow  a  rapid  relaxation  by 
losing  energy  to  zone  edge  modes  through  the  anharmonic  forces.  One 
interesting  result  of  the  ultrafast  relaxation  experiments  has  been 
the  unexpectedly  slow  cooling  of  the  hot  distribution.  The  cooling 
rate  has  generally  been  found  to  decrease  with  increasing  carrier 
densities.  In  the  case  of  a  far-from-equilibrium  electron-hole 
plasma,  one  would  expect  to  see  the  consequences  of  production  of 
excess  phonons.  The  hot  phonon  effect  emphasizes  the  excess  phonon 
population,  while  the  essential  issue  of  the  problem  is  phonon  reab¬ 
sorption  by  the  carriers  which  really  slows  the  cooling  rate.  The 
probability  of  phonon  reabsorption  increases  with  increasing  carrier 
density.  The  lifetime  of  LO  phonons  has  been  measured  to  be  ap¬ 
proximately  7  psec  in  GaAs  [6~8).  If  the  carrier  density  is  large 
enough,  the  optical  phonon  will  be  reabsorbed  before  it  has  chances  to 
decay.  Ve  assume  that  once  the  phonon  decays,  the  energy  is  lost. 
Von  der  Linde  et_  al_.  [ 6 j  showed  that  the  anharmonic  decay  of  the  LO 
phonons  essentially  explains  the  lifetime  in  bulk  GaAs. 

Free  carrier  screening  has  also  been  confirmed  to  reduce  the 
cooling  rates  to  some  extent.  Nevertheless,  it  still  is  not  suffi¬ 
cient  to  interpret  the  experimental  results.  However,  including  hot 
phonon  effects  has  led  to  a  much  closer  agreement  [19].  Screening 
will  affect  which  modes  become  hot,  and  the  degree  to  which  the  phonon 
population  is  disturbed.  Furthermore,  screening  is  very  effective  at 
reducing  the  scattering  of  electrons  by  the  Frohlich  interaction  since 
this  interaction  is  strongest  at  small  q-vectors,  where  screening  is 


most  significant.  Screening  has  less  influence  on  holes  than 


Since  the  initial  energy  of  the  electrons  is  much  higher  than  the 


holes,  the  electron-hole  plasma  forms  a  f ar-f rom-equilibrium  system. 
One  would  expect  the  energy  to  be  transferred  from  the  hot  electrons 
to  the  cold  holes  and  to  the  lattice,  through  the  electron-hole  inter¬ 
action  and  carrier-phonon  interaction,  respectively.  At  low 
excitation  levels  the  energy  flow  is  dominated  by  the  electron-phonon 
interaction,  while  at  high  excitation  levels  the  energy  is  mainly 
channelled  via  electron-hole  interaction  to  the  holes  and  then  via  the 
TO  phonon  emission  by  the  holes,  due  to  the  unscreened  deformation 
potential  interaction  with  the  lattice.  This  is  understandable, 
because  at  high  carrier  densities  the  electron-phonon  interaction  is 
highly  screened  while  the  electron-hole  interaction  becomes  stronger. 
Since  the  electron-hole  scattering  rate  is  high,  a  large  fraction  of 
each  electron's  energy  will  eventually  be  channelled  to  the  holes  when 
an  equilibrium  is  reached  between  the  electron  and  hole  systems.  In 
addition,  the  electron-hole  interaction  slows  the  cooling  rates  at  low 
carrier  densities  by  gradually  shifting  the  electron  population  to  the 
lower  energy  states.  However,  it  enhances  the  cooling  rate  at  high 
carrier  densities  because  the  electrons  can  transfer  energy  to  the 
lattice  through  the  holes.  Fig.  2.10  fully  depicts  the  energy  flowing 
channels  for  laser-excited  system  at  different  carrier  densities.  A 
first  experiment  investigating  the  electron-hole  energy  transfer  in 
semiconductors  has  been  reported  by  Hbpfel  and  coworkers  [14] 
recently.  They  observed  that  the  energy-loss-rate  of  hot  minority 


electrons  in  a  cold  hole  plasma  is  higher  by  about  a  factor  of  two 
than  that  of  majority  electrons  in  the  absence  of  holes. 


CHAPTER  3 


ENSEMBLE  MONTE  CARLO  TECHNIQUES  FOR 
HOT  PHOTOEXCITED  CARRIER  RELAXATION 

3.1  Boltzmann  Transport  Theory 

As  a  preliminary  to  any  statistical  discussion  of  nonequilibrium 
processes,  one  must  find  a  vay  to  evaluate  the  distribution  function 
when  the  system  is  nonequilibrium.  This  is  accomplished  by  means  of 
the  Boltzmann  equation.  If  the  rate  of  change  of  the  distribution 
function  due  to  collisions  or  scatterings  is  denoted  by  c 0 1 1 ’  t^ien 
the  total  rate  of  change  of  the  distribution  function  f  can  be 
expressed  as 


df_  _  3f  _  ej  3f  ran 

dt  =  dr  '  h  'dit  +  ^co11’ 


(3.1) 


where  f  is  the  applied  external  field  and  mv  =  ftk  is  the  momentum. 
The  first  and  second  term  on  the  right  hand  side  of  (3.1)  represent 
the  effect  of  the  density  gradient  and  of  the  applied  external  field 
on  the  distribution  function.  In  semiconductors ,  the  behavior  of 
carriers  is  affected  by  the  applied  external  field,  density  gradient, 
and  scattering  from  phonons,  impurities  or  other  carriers  etc.  All 
these  factors  can  be  included  in  the  Boltzmann  transport  equation  to 
obtain  information  about  the  carrier  distribution  function  f.  This 
semi -classi cal  integro-di f ferential  equation  for  f  can  be  solved 


numerically  with  the  help  of  fast  computers.  Otherwise,  some 
simplifying  assumptions  have  to  be  made  in  order  to  solve  the  equation 
analytically  under  certain  special  conditions.  The  most  generally 
used  analytical  methods  are  the  Legendre  polynomials  expansion  and  the 
displaced  Maxwellian  approximation  [47].  One  of  the  numerical 
techniques  is  the  iterative  method,  while  the  other  one  is  the  Monte 
Carlo  method  on  which  we  focus  our  attention. 

3.2  Application  of  Ensemble  Monte  Carlo  Method  to  Transport  in 

Semiconductors 

The  Monte  Carlo  method  was  first  applied  to  the  transport  in 
semiconductors  by  Kurosawa  (1966)  to  study  steady-state  hole  transport 
in  Ge.  It  is  a  very  powerful  numerical  technique,  because  it  can 
simulate  special  physical  situations  which  are  unachievable  in 
experiments,  and  also  investigate  the  non-existing  materials  to 
examine  important  features  of  phenomenon  in  vthich  we  are  interested. 
Since  this  stochastic  evaluation  of  the  solution  of  the  Boltzmann 
transport  equation  is  a  direct  way  to  simulate  the  dynamics  of  the 
carriers  inside  semiconductors,  any  physical  information  can  be  easily 
obtained  while  the  time-evolving  solution  of  the  equation  is  being 
determined.  The  single  particle  Monte  Carlo  method  is  extended  to  a 
many-particle  ensemble  Monte  Carlo  method  in  order  to  investigate  non¬ 
steady-state  or  non-homogeneous  phenomena  in  semiconductors.  A 
general  flow  chart  of  the  ensemble  Monte  Carlo  program  is  given  in 


The  semiclassical  EMC  method  contains  simulations  of  the  motion 


of  an  ensemble  of  charge  carriers  inside  the  semiconductor  vith  (or 

without)  an  electric  field  applied  and  vith  the  scattering  mechanisms 

specified.  The  scattering  probabilities  of  each  scattering  mechanism 

are  calculated  and  tabulated  inside  subroutines  of  the  program.  The 

alternative  sequence  of  drifts  and  instantaneous  collisions  vith 

phonons,  impurities  or  other  carriers  are  iteratively  simulated, 

relying  on  the  generation  of  a  sequence  of  random  numbers.  A  special 

method  to  determine  the  free  flight  time  without  solving  a  complicated 

integral  equation  directly  vas  devised  by  Rees  [54],  vho  introduced  a 

fictitious  "self-scattering"  so  that  the  total  scattering  probability, 

including  this  self  scattering,  is  equal  to  a  constant  T  which  is 

taken  to  be  at  least  as  large  as  the  largest  value  of  the  total 

scattering  rates.  Then,  P0(E)  =  T  —  P(E)  is  the  probability  of  self 

scattering.  If  a  carrier  undergoes  such  a  self  scattering  which  is  of 

no  physical  significance  and  does  not  change  the  distribution 

— * 

function,  its  final  state  after  collision  is  the  same  as  the  initial 
state  before  the  scattering,  so  that  the  electron  motion  remains 
unperturbed  as  if  no  scattering  occurred  at  all.  A  random  number  r 
distributed  between  0  and  1  can  be  used  to  generate  a  stochastic  free 
flight  time  t ^ ,  vhere 


i 


It  is  advantageous  to  make  T"  as  small  as  possible  in  order  to  minimize 
the  number  of  self-scattering  events  [  9  J  .  The  computer  time  wasted  ir. 


taking  care  of  self-scat tering  is  more  than  compensated  by  the 
simplification  of  the  calculation  required  to  determine  the  free 
flight  time.  Having  determined  the  free  flight  time,  it  is  necessary 
to  determine  which  scattering  process  terminates  the'  free  flight.  The 
arbitrary  choice  of  the  scattering  mechanism  for  a  particular  carrier 
with  energy  E  is  also  determined  by  the  random  numbers.  Since  the 
probability  of  the  electron  being  scattered  by  process  i  is 
proportional  to  the  scattering  rate  P^(E)  for  the  ith  mechanism,  and 
since  ^qP^(E)  =  T,  it  is  necessary  to  generate  a  random  number  r 

between  0  and  1  to  test  the  inequality  rT  <  .5^P^(E)  for  all  m.  When 
this  inequality  is  satisfied  the  mth  scattering  process  is  chosen. 
Figure  3.2  illustrates  how  a  mechanism  is  picked,  as  elucidated.  If 
all  scattering  processes  have  been  tried  and  none  has  been  selected, 
it  means  that  rT  >  P(E),  and  a  self  scattering  occurs,  where  P(E)  = 
.^P.fE).  The  energy  of  the  final  state  is  determined  by  the 
properties  of  the  chosen  scattering  process,  while  the  determination 
of  the  orientation  of  the  final  wave-vector  requires  further 
generation  of  random  numbers  to  determine  the  angles  G  and  4>  (see  Fig. 
3.3)  according  to  the  angular  dependence  of  the  chosen  scattering 
process  described  in  Chapter  2. 

This  EMC  method  traces  the  history  of  an  ensemble  of  carriers,  in 
parallel,  for  successive  very  short  time  intervals,  instead  of 
following  the  path  of  a  single  particle  for  a  long  period  of  time. 


Each  carrier  in  the  ensemble  can  be  described  by  a  set  of  physical 
variables  (Q.}  that  fully  describe  the  carrier's  state.  Since  every 


Q.  is  calculated  at  each  time  step  using  the  MC  method,  the  set  {0.} 
is  also  an  ensemble  evolving  in  time.  The  average  value  of  a  quantity 
Q  is  defined  as  the  ensemble  average  at  time  t  over  the  N  carriers  of 
the  system, 


<Q>  =  Nf0i(t>* 


3.3  Ensemble  Monte  Carlo  Model  for  Electron-Hole  Interaction 


It  is  very  hard  to  evaluate  the  carrier-carrier  scattering  rates 

analytically  because  one  has  to  assume  the  shapes  of  the  distribution 

functions  f  ,  f,  involved  in  the  integrals  as  expressed  in  (2.4)  and 
e  h 

(2.5).  However,  one  of  the  beauties  of  the  EMC  technique  is  that, 
there  is  no  need  to  manipulate  the  complex  integrals  analytically 
since  the  scattering  rates  can  be  expressed  as  the  ensemble  averages 
in  the  following, 


Peh(kl>  = 


Fhe<kl>  = 


oue  l 


2  nc^h^  Nh1  g2  (o^h-6"  j 
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where  K, ,  N  are  the  total  numbers  of  electrons  and  holes  in  the 
h  e 


ensemble  respect ivelv ■ 


3.3.1  Description  of  Flow  Chart  of  the  Model 


A  flow  chart  of  the  EMC  program  developed  to  study  the  relaxation 
of  the  hot  photoexcited  carriers  is  illustrated. in  Fig.  3. A.  The 
initialization  step  gives  the  material  parameters,  the  initial 
distribution  of  carriers,  and  other  input  parameters  needed  to 
investigate  the  relaxation  such  as  the  number  of  iterations,  the 
observing  time,  lattice  temperature,  the  photon  energy  and  pulse 
duration,  the  flags  to  specify  particular  conditions  of  interest,  etc. 
The  second  step  is  to  calculate  and  tabulate  the  scattering 
probabilities  for  each  mechanism,  such  as  acoustic,  polar  optical, 
non-polar  optical  scatterings,  etc.  which  are  all  calculated  without 
screening  in  this  step.  The  third  step  is  photoexcitation.  During 
this  step,  the  electron-hole  pairs  are  generated  according  to  the 
given  pulse  shape  and  the  energies  are  distributed  to  the  carriers 
according  the  band  structure.  Also,  the  spatial  distributions  of  the 
photo-generated  carriers  are  determined.  The  fourth  block  is  to 
calculate  the  carrier-carrier  scattering  rate  using  the  formulae 
developed  above.  Simultaneously,  tr.e  sei  f -cons  i  s  t  er.  t  screening  length 
is  also  determined  for  calculating  the  screened  carrier-phonor. 
scattering  rates.  The  fifth  block  :r  • ?  : a  1  - c 1  a  *  e  'he  screened 
carrier-phonon  scattering  rater  using  .ng  1  g::.  caicula'er 
previously.  The  dynamics  •:  :  •  •  •  : «  .v.-f  *  i gates, 
using  EMC  method  sepe:  a*<  lv. 

EMC.  because,  once  the  '  ■  .  •  . . 
be  selected  from  the  >•••■■•■  ■  •  «  :  k- c 


coupled  to  the  electron  EMC  because  of  the  hole-electron  scattering. 
It  is  noted  that  the  scattering  rates  of  electron-hole  scattering  and 
hole-electron  scattering  are  different.  In  the  eighth  block,  the 
ensemble  averages  of  interest  are  calculated,  such  as  the  mean  energy 
of  carriers,  the  equivalent  carrier  temperature,  or  drift  velocity 
etc.  Additionally,  the  energy  and  momentum  distributions  of  carriers 
are  specified  in  this  block  too.  When  the  pulse  is  still  on,  steps 
3~8  are  repeated.  Only  steps  4-8  are  repeated  when  the  pulse  is  over. 

3.3.2  Self-Consistent  Screening  Model 

Apparently,  the  screening  length  has  great  influence  on  the 
screened  carrier-phonon  and  carrier-carrier  scattering  rates.  Several 
analytical  models  for  calculating  the  screening  length  have  been 
reported  shoving  that  the  magnitude  of  the  screening  length  varies 
vith  different  models.  A  self-consistent  screening  length  in  the 
static  and  long  vavelength  limit  of  the  random  phase  approximation 
calculated  in  the  EMC  method,  taking  advantage  of  its  built-in  time- 
evolving  distribution  functions,  has  been  developed  recently  [35]. 
The  vay  it  is  calculated  is  described  in  the  following.  The 
dielectric  constant  e(q,w)  is  given  by 

,  ,  .  e2  „  f (k) 

C(q’w)  '  1  2  £  E(k*q) 

t  a 

CD  * 

In  the  static  and  long  vavelength  case,  (3.2)  becomes  [56] 
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e(q,0)  -  1  ♦  - -l\-  — ^-jD(E)dE 

C«fl 


where  D(E)  is  the  density  of  states  and  6  is  the  inverse  of  screening 
length  given  as 


af(E)') 

“SE  JD(E)dE 


(3.3) 


3D(E) 

3E 


f (E)dE. 


(3.4) 


For  convenience,  ve  define  B(E)  as 


=  B(E)D(E) . 


(3.5) 


Obviously,  B(E)  is  the  logarithmic  derivative  of  the  density  of  states 
D(E).  Inserting  (2.5)  into  (3.4),  ve  obtain 


r  =  §-  J  B(E)D(E)f (E)dE, 

OD 


B(E)> 


(2.6) 


where 


ill 


<B(E)>  = 


B(E)D(E)f(E)dE 
J  D(E)£(E)dE  ' 


and  the  carrier  concentration  n  is  defined  as 


n  =  /  D(E) f (E)dE. 


From  (3.6),  ve  see  an  amazing  result  shoving  that  g2  can  be  expressed 
in  terms  of  the  ensemble  average  of  B(E)  and  is  given  as 


g2  =  —  J-  r  b(e . ) , 

H  c  N  i  1 


(3.7) 


where  N  is  the  number  of  carriers  in  the  ensemble.  The  form  of  B(E) 
depends  on  the  band  model  used  in  the  EMC  calculation.  For  parabolic 
bands,  B(E)  is  given  as 


B(E)  = 


For  nonparabolic  bands,  B(E)  becomes 


B(E)  = 


l-'-ZaE 


l*2c£ 

2E(l*o£) 


Vhile  assuming  a  Kaxvell-Boltzmann  distribution  function,  (3.6) 


becomes 
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which  is  the  screening  length  of  the  Debye-HUckel  screening  model. 
and  are  the  electron  and  hole  temperatures,  respectively.  The 

cooling  rate  predicted  by  the  Debye-HUckel  model  is  slower  than  the 
one  calculated  by  the  self-consistent  screening  model,  which  means  the 
Debye-HUckel  model  overestimates  the  screening.  The  results  of  the 
comparisons  will  be  discussed  in  the  next  chapter. 

3.4  Ensemble  Monte  Carlo  Model  for  Hot  Phonon  Effects 

In  addition  to  the  difficulties  needed  to  be  overcome  in 
programming,  the  electron-electron  and  hole-hole  interactions,  and  the 
screening  of  carrier-phonon  interactions,  as  well  as  electron-plasmon 
scatterings  are  all  ignored  in  this  simulation  to  better  demonstrate 
the  hot  phonon  effects  and  electron-hole  interactions.  The  time- 
evolving  LO-phonon  distribution  function  is  specified  and  the  modified 
scattering  rates  of  the  electrons  due  to  the  nonequilibrium  phonon 
population  are  recalculated.  Ve  only  consider  the  disturbance  of  the 
LO-phonons  while  ignoring  the  other  types,  because  the  LO-phonons  play 
a  major  role  in  most  of  the  cooling  process.  The  algorithm  developed 
to  study  the  hot  phonon  effects,  which  slows  down  the  electron 
relaxation  in  the  bulk  laser-excited  GaAs,  will  be  discussed  below. 

First,  we  need  to  set  up  a  histogram  h^  defined  over  the  q  space 
for  the  LO-phonons  in  the  program  to  follow  the  time-evolving  phonon 
distribution.  After  each  time  step,  in  which  the  electron  suffers  a 
collision  involving  an  LO-p'nonon,  the  histogram  is  updated.  At  fixed 
successive  time  intervals  ( n At ) ,  where  At  is  always  much  shorter  than 


I 
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the  average  scattering  time  for  the  LO-phonon  scattering,  can  be 
calculated  by  [25] 


N  (nAt)  =  N  (nAt)  -  |N  (n  At  )-N0 1-^— , 
q  q  L  q  J  T 


(3.8) 


vhere 


Nq(nAt)  =  l (n-1) At ]+AAh^ , 


and  N0  is  the  equilibrium  Bose-Einstein  phonon  distribution  function. 
The  term  AAh^  gives  the  dynamical  contribution  of  the  electronic 
processes  to  the  phonon  distribution  during  the  time  step  At,  vhere  A 
is  a  normalization  factor  accounting  for  the  density  of  states  in  q 
space  and  for  the  concentration  of  excited  electrons.  A  is  given  as 


Nq  ( Aq) 


vhere  n  is  the  concentration  of  electrons  in  the  central  valley,  N  is 
the  total  number  of  electrons  in  the  ensemble,  q  is  the  phonon  vave 
vector,  and  Aq  is  the  incremental  grid  element  of  phonon  vave  vector 
in  the  histogram.  The  second  term  on  the  right  hand  side  of  (3.8) 
represents  the  phonon-phonon  processes  vith  an  experimentally 
determined  LO-phonon  relaxation  time  ~  of  7  ps  16],  In  Lugli's  paper 
[25]  ,  the  scattering  rates  concerning  the  LO-phonons  are  recalculated 
at  every  fixed  time  interval  (nAt)  according  to  the  nev  time-evolving 
phonon  distribution  function  given  by  (3.8),  so  that,  the 
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modifications  resulting  from  the  phonon  disturbance  are  taken  into 

account.  Here,  instead  of  recalculating  the  scattering  rates  for  the 

LO-phonons  using  the  perturbed  phonon  distribution,  a  very  large 

number  chosen  at  the  beginning  of  the  simulation  has- been  used  for  the 

maximized  nonequilibrium  phonon  distribution  function  to  calculate  the 

scattering  rates  for  LO-phonons  and  then  a  rejection  method  is  used  to 

determine  the  final  state  after  each  scattering  process  involving  an 

LO  phonon  [57].  The  algorithm  expressed  above  has  been  used  to  study 

0 

the  relaxation  process  of  the  bulk  GaAs  at  nitrogen  temperature  (77  K) 
illuminated  by  a  laser  pulse  of  either  delta  pulse  shape  or  0.2  ys 
pulse  duration.  The  results  of  the  calculations  will  be  discussed  in 
next  chapter. 


CHAPTER  4 


I 

i  RESULTS  AND  DISCUSSIONS 

j  The  influences  of  the  electron-hole  interaction  and  hot  phonon 

effects  to  the  ultrafast  relaxation  of  hot  laser-excited  carriers  in 

'  GaAs  vere  investigated  via  the  ensemble  Monte  Carlo  method  for 

16  18  -3 

|  different  excitation  levels,  ranging  from  5x10  to  1x10  cm  ,  at 

different  excitation  energies.  The  pulse  duration  vas  taken  to  be  a 
delta  function  for  convenience  in  the  calculations.  Additionally,  the 
cooling  rates  of  the  photoexcited  carriers,  while  disregarding 
carrier-carrier  interactions,  were  evaluated  using  both  the  Debye- 
HUckel  and  the  self-consistent  screening  models  which  have  already 
been  built  in  the  program.  The  results  of  the  simulations  will  be 
discussed  and  compared  in  the  following  sections. 

4.1  The  Effect  of  Electron-Hole  Interaction 

The  initial  mean  energy  of  the  excited  electrons  is  about  seven 
times  the  corresponding  energy  of  holes,  according  to  the  band 
structure  of  GaAs,  so  that  the  electron-hole  plasma  forms  a  far-from 
equilibrium  system.  Accordingly,  the  energy  is  transferred  from  the 
hot  electrons  to  the  cold  holes  and  then  to  the  lattice  through  the 
electron-hole  and  carrier-phonon  scatterings.  In  this  section,  the 
manner  in  which  the  electron-hole  interaction  affects  the  role  of  the 
electron-phonon  interaction.  which  dominates  the  cooling  processes, 
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was  examined.  The  phonons  were  assumed  to  be  in  thermal  equilibrium 
in  this  situation. 

The  time  evolution  of  the  electron  mean  energy  is  shown  in  Fig. 
4.1  to  illustrate  hov  the  electron-hole  interaction  behaves  under 
different  carrier  concentrations.  It  is  found  that  the  cooling  rate 
slows  down  as  the  excitation  level  increases  under  the  situation  where 
only  unscreened  carrier-phonon  and  electron-hole  interactions  were 
active.  The  holes  are  expected  to  behave  the  same  way.  The  relative 
importance  of  the  electron-phonon  and  electron-hole  interactions  as 
energy  transfer  channels  depends  on  the  excitation  level  and 
excitation  energy.  From  Fig. 4. 2,  it  is  found  that  the  electron-hole 
interaction  promotes  the  cooling  rate  at  high  excitation  levels 
because  the  hot  electrons  transfer  their  excess  energy  to  the  lattice 
through  the  holes  mainly  by  emission  of  TO  phonons  via  the  unscreened 
optical  deformation  potential  coupling.  However,  the  electron-hole 
interaction  slows  down  the  cooling  rate  at  low  excitation  level  as 
shown  in  Fig.  4.3  because  the  electron  population  was  shifted 
gradually  to  lower  energy  states  via  the  electron-hole 
scat  terings [ 56 ] . 

It  is  obvious  that  the  electrons  cool  faster  at  lover  excitation 
levels  as  shown  in  Fig.  4.1  where  the  screening  is  negligible.  The 
slowing  becomes  more  significant  at  high  carrier  densities,  where  the 
optical  phonons  are  strongly  screened  and  the  electron-electron 
scatterings  occur  more  frequently.  The  electron-electron  and  hole- 
hole  interactions  were  ignored  in  this  calculation  because  of  the 
complications  of  some  major  changes  in  the  program.  However,  it  is 
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anticipated  that  the  hole-hole  scatterings  enhance  the  cooling  of  the 
hole  system  while  the  electron-electron  scatterings  slow  the  cooling 
of  the  electrons.  The  reason  why  the  cooling  of  the  holes  is  enhanced 
can  be  resorted  to  the  fact  that  the  hole-hole  scattering  rate  is  very 
high  so  that  the  hole-hole  scatterings  kick  some  holes  above  the 
phonon  emission  threshold,  and  these  then  emit  phonons  through  the 
deformation  potential  interaction.  On  the  other  hand,  the  interaction 
among  the  electrons  quickly  redistributes  the  electrons  into  high  and 
lov  energy  regions[36]. 

4.2  The  Role  of  Hot  Phonon  Effects 

The  model  used  to  examine  the  hot  phonon  effects  is  basically 
based  on  the  model  proposed  by  Lugli  and  coworkers [ 25 ] .  This  model 
takes  account  of  the  electron-hole  and  non-equilibrium  LO  phonon- 
electron  interactions.  Electron-electron,  hole-hole,  and  electron- 
plasmon  interactions  and  the  screening  of  carrier-phonon  interactions 
were  ignored  here.  The  effect  of  hot  phonons  was  examined  for  an 
excitation  energy  of  1.8  ev,  which  does  not  exceed  the  band  gap  of 
1.51  eV  too  much,  in  order  to  prevent  complications  caused  by  exciting 
electrons  to  the  upper  valleys.  It  is  found  that  the  non-equilibrium 
phonons  do  slow  down  the  cooling  rate  of  the  hot  laser-excited 
electrons  further,  as  shown  in  Fig.  4.4,  by  comparing  with  the  results 
obtained  from  the  situation  where  only  the  unscreened  carrier-phonon 
and  the  electron-hole  interactions  were  active.  This  phenomenon  is 
mainly  caused  by  the  reabsorption  of  previously  emitted  LO  phonons 


63 


without  having  enough  time  to  decay  into  near-zone  boundary  phonons, 
which  would  not  couple  back  to  the  electrons,  as  well  as  the 
thermalization  of  the  electrons  and  holes.  In  Fig.  4,4,  about  1  ps 
after  the  laser  pulse,  phonon  heating  effects  become  dominant  for 
reducing  the  cooling  rates. 

The  time  evolving  LO  phonon  distribution  function  N(q)  is 

calculated  in  terms  of  the  phonon  wavevector  q  and  assumed  to  be 

isotropic  in  q.  Figures  4.5-8  trace  the  non-equilibrium  LO  phonon 

0 

populations  evolving  with  time.  At  77  K,  the  value  of  the  equilibrium 

-3 

phonon  distribution  is  about  4.179x10  ,  as  indicated  by  the  arrow  in 

Fig.  4.5.  It  is  known  that  the  strongly  perturbed  phonon  distribution 
is  due  solely  to  the  strong  LO-phonon  emission  of  the  hot  electrons. 
After  the  excitation,  reabsorption  of  the  emitted  phonons  by  the 
carriers  and  the  phonon-phonon  losses  bring  the  distribution  back  to 
the  equilibrium  value.  It  is  noted  that  even  a  few  picoseconds  after 
the  laser  pulse,  a  significant  number  of  phonons  still  exist,  and  a 
considerable  number  of  phonon  reabsorptions  still  go  on.  Since  the 
numbers  of  electrons  and  holes  residing  in  the  ensemble  being 
simulated  are  taken  to  be  as  small  as  2000  respectively  in  order  to 
save  the  computer  time,  the  crude  approximation  results  in  the  jumps 
of  the  phonon  distribution  functions,  as  shown  in  Figs.  4.5-8,  which 
are  expected  to  be  smooth  curves.  The  disturbances  of  other  types  of 
phonons  are  ignored  in  this  calculation.  Actually,  emission  of  TO 
phonons,  which  only  couple  to  the  holes,  also  provides  a  contribution 
to  the  energy  loss [16]. 
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As  expected,  hot  phonon  effects  are  mainly  responsible  for  the 
relatively  high  carrier  temperatures  observed  in  experiments 
investigated  in  the  subpicosecond  time  regime.  The  deviation  of  the 
phonon  population  from  its  equilibrium  value  further  slows  down  the 
hot  electron  relaxation,  which  bring  closer  agreement  between 
theoretical  and  experimental  results  for  laser-excited  GaAs . 

4.3  Screening  Effect  of  Electron-Phonon  Interaction 

Free  carrier  screening  reduces  the  cooling  rates,  but  not 
sufficiently  to  explain  the  experimental  data.  It  is  noted  that 
inclusion  of  hot  phonon  effects  leads  to  closer  agreement  between 
experiment  and  theory.  In  this  section,  the  cooling  rates  of  the 
laser-excited  electrons  were  calculated  in  the  absence  of  the  carrier- 
carrier  interactions,  hot  phonon  effects  and  in  the  presence  of  holes 
using  both  the  Dcbye-HOckel  and  the  self-consistent  screening  models 
to  better  demonstrate  the  screening  effcet. 

Figure  4.9  illustrates  that  the  screening  effect  reduces  the 
cooling  rates,  as  compared  to  the  case  where  screening  is  ignored. 
The  lowest  curve  represents  the  cooling  rate  when  screening  is  totally 
ignored,  which  is  too  fast  compared  to  the  observed  slow  cooling  rates 
at  high  excitation  levels.  Vhen  the  screening  of  the  electron-phonon 
interaction  is  considered,  the  cooling  rate  is  reduced  to  a  different 
extent  through  the  two  different  screening  models.  In  the  self- 
consistent  screening  model.  the  screening  evolves  as  the  carriers 
relax  and  the  screening  length  is  calculated  self-consistently  as 
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discussed  in  Sec.  3.3.2.  The  reduction  of  the  cooling  rate  is  smaller 
for  the  self-consistent  model,  which  gives  more  reasonable  estimates 
of  the  screening  length  and  reflects  the  time  evolution  of  the  carrier 
distribution  functions.  The  Debye-HUckel  model  overestimates  the 
screening  due  to  the  equilibrium  Boltzmann  distribution  function 
imposed  on  the  carriers. 

4.4  Time  Evolution  of  Electron  Energy  Distribution  Function 

In  this  section,  the  manner  in  which  the  electron-hole  and 

carrier-phonon  interactions  affect  the  form  of  the  energy  distribution 

function  of  hot  electrons  is  examined.  Figure  4.10  illustrates  the 

distribution  function  of  electrons  at  0.25  ps  after  a  1.8  ev,  delta- 
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shaped  laser  pulse  at  an  excitation  level  of  5x10  cm  .  This  figure 
displays  peaks  below  the  initial  energy,  displaced  at  integral 
multiples  of  the  L0  phonon  emission  threshold  (about  36.4  mev).  Even 
at  1  ps  after  the  laser  pulse,  the  distribution  function  still  shows 
some  peaks,  as  depicted  in  Fig.  4.11,  which  assures  us  that  the 
electrons  cascade  down  to  lover  energies  mainly  by  emitting  LO 
phonons.  It  is  shown  that  the  electron-phonon  interaction  dominates 
the  cooling  at  low  excitation  levels.  The  peaks  are  broadened  by  the 
electron-hole  scatterings  and  it  seems  that  there  are  more  electrons 
on  the  lower  side  of  the  peaks,  indicating  the  energy  loss  by  the 
electrons  during  the  electron-hole  scatterings.  Figure  4.12 
illustrates  the  distribution  function  at  5  ps,  which  seems  becoming  a 
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CHAPTER  5 


CONCLUSIONS 

The  ensemble  Monte  Carlo  technique  provides  very  accurate 
microscopic  information  about  the  physical  system  under  simulation  and 
it  does  not  require  any  assumption  on  the  carrier  distribution 
function  because  of  the  self-consistent  time-evolving  distribution 
function  inherent  vithin  it.  Additionally,  using  a  more  realistic 
band  structure  and  a  direct  evaluation  of  the  scattering  probabilities 
for  the  various  scattering  mechanisms  involved  in  the  relaxation 
process  becomes  attainable. 

The  major  purpose  of  this  thesis  is  to  develop  an  ensemble  Monte 
Carlo  method  to  simulate  the  hot  laser-excited  electron-hole  system  in 
bulk  GaAs,  taking  account  of  the  carrier-phonon  interaction,  electron- 
nole  interaction,  and  hot  phonon  effects  in  order  to  obtain  closer 
agreement  betveen  the  theoretically  calculated  and  experimentally 
observed  cooling  rates  of  hot  p'notoexcited  carriers  in  semiconductors. 
Generally,  the  relaxation  of  hot  laser-excitec  carriers  can  be  divided 
into  three  stages.  The  first  stage  is  the  \ery  fast  initial 
relaxation  vhere  the  carriers  lose  most  of  their  excess  energy.  In 
the  second  stage,  the  hot  distribution  cools  tovards  the  lattice 
temperature  until  the  optical  phonons  are  no  longer  able  to  remove 
energy.  The  first  tvo  stages  take  typically  over  vithin  100  ps.  The 
third  stage  occurs  in  the  nanosecond  tire  region  involving  the  final 
cooling  of  the  distribution  to  the  lattice  temperature  through  the 
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It  is  known  that  the  electron-phonon  interaction  is  the  dominant 
channel  for  energy  loss  by  hot  electrons  at  low  excitation  levels, 
while  at  high  excitation  levels,  the  electron-hole  interaction  is  the 
main  channel.  The  holes  then  dissipate  energy  through  the  hole-phonon 
interaction  via  the  unscreened  TO  phonon  deformation  potential 
interaction.  The  electron-hole  interaction  promotes  the  cooling  rate 
of  hot  electrons  at  high  carrier  densities  by  channeling  the  electron 
energy  to  the  holes  and  then  to  the  lattice.  However,  it  slows  down 
the  cooling  rate  at  low  carrier  densities  due  to  the  shifting  of  the 
electron  population  to  energies  below  the  LO  phonon  emission 
threshold.  The  cooling  rate  has  been  found  to  decrease  with 
increasing  excitation  level  because  both  hot  phonon  effects  and 
screening  become  more  important  in  a  high  density  plasma.  If  the 
carrier  density  is  large  enough,  the  optical  phonons  will  be 
reabsorbed  before  they  start  to  decay.  The  experimentally  observed 
high  kinetic  energy  of  the  electron-hole  plasma  is  sh  :  -’n  to  be  mainly 
due  to  LO  phonon  heating  which  provides  the  dominant  mechanism  to 
explain  the  slowing  down  of  the  cooling  rate  as  well  as  the 
thermalization  of  the  electrons  and  holes. 

Recently,  there  has  been  an  extensively  interest  in  the  two 
dimensional  quantum  well  structure.  The  behavior  of  the  electrons  in 
quantum  veil  differ  drastically  from  the  three  dimensional  bulk 
semiconductors.  A  further  understanding  of  carrier  cooling  in  quantum 
well  structures  requires  more  detailed  knowledge  of  the  phonons  n 
these  systems.  Extending  the  present  work  to  two  cimensional  quantum, 
well  structures  will  be  an  exciting  task  to  carry  on  in  the  near 
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